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We present a theoretical study of electron-phonon scattering effects in thin films made of a strong 
topological insulator. Phonons are modelled by isotropic elastic continuum theory with stress-free 
boundary conditions, and the interaction with the helical surface Dirac fermions is mediated by the 
deformation potential. We determine the temperature-dependent electrical resistivity p(T) and the 
quasiparticle decay rate F(T) observable in photoemission. The low- and high-temperature power 
laws for both quantities are obtained analytically. Detailed estimates covering the full temperature 
range are provided for Bi2Se3. 
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I. INTRODUCTION 

The recently discovered state of matter called "topo- 
logical insulator" (TI) currently represents one of the 
most active areas in condensed matter physics* 1 ^ TIs are 
characterized by an insulating gap in the bulk but at the 
same time have an odd number of gapless surface modes 
protected against all time-reversal invariant (and suffi- 
ciently weak) perturbations?^ 4 - In a three-dimensional 
(3D) TI, these surface modes correspond to massless two- 
dimensional (2D) Dirac fermions, where the spin direc- 
tion is in the surface plane and perpendicular to momen- 
tum ("spin-momentum locking"). A typical reference 
material is Bi2Se3 with a bulk gap Af, ss 0.3 eV. The 
helical Dirac electron property of the TI surface state 
has been experimentally confirmed by spin- and angle- 
resolved photoemission spectroscopy (ARPES).-— Trans- 
port experiments arc more difficult in that respect since 
the surface contribution is often masked by the residual 
conductivity due to impurities or defects in the bulk£~— 
In thin films made of TI materials, however, the bulk 
contribution is largely suppressed relative to the surface 
contribution, rendering the latter easier to observe. 

In this paper, we provide a detailed theoretical 
analysis of both the temperature-dependent resistivity 
p(T) and the quasiparticle lifetime T(T) (observable in 
ARPES2ii£) for a thin TI film. The approach taken here 
generalizes previous work for the semi-infinite geometry 
(with only one surface) by two of usiL to the film geom- 
etry. This brings about several important changes com- 
pared to Ref. [ll| that are discussed below. We model the 
electronic part by retaining only the Dirac surface states 
obtained from the low-energy bandstructure^ and our 
theory always assumes that the Fermi level is located 
inside the bulk gap. We note in passing that Ref. UH pro- 
vides more accurate parameter estimates as the earlier 
paper by the same authors^ and we here adopt their 
new parameters in our calculations using Bi2Se3 as ex- 
ample. A similar parameter set has been published in 
Ref. OJ. 

In sufficiently thin films, the hybridization of the two 
surface states eventually causes insulating behavior, as 
has recently been observed experimentally from ARPES 



for Bi2Se3 films . 15 ' 16 For Bi2Se3, several calculations 
predicte d 14 ' 17 ' 18 a gap A(i) with (as a function of the 
width L) oscillations superimposed on an exponential 
decay. Similar calculations, however, found no oscilla- 
tions, with a well-established TI phase already for L > 3 
quintuple layers (QLs)i 19 ' 20 Using the parameter set of 
Ref. KL21 we also find no evidence for oscillations in A(L), 
see Sec. Ill Al below. For large width, one then has (upper 
and lower) massless Dirac fermion surfaces? 2 ^ 

Our working assumption below is that electron-phonon 
scattering is the dominant source of quasiparticle de- 
cay and backscattering. Electron-electron interactions 
are indeed expected to give only subleading corrections 
to the resistivity as long as T > 1 niK? 2 * 2 . Disorder ef- 
fects are more likely to compete with phonon-induccd 
backscattering effects. However, for elevated tempera- 
tures, T > 100 K, phonon effects dominate even for 
present-day samples, and anticipating higher purity films 
in the future, this crossover temperature may be low- 
ered significantly. ARPES setups allowing to test our 
predictions for the quasiparticle decay rate are basically 
available ££2r2k Other surface scattering techniques have 
also been applied to extract the phonon dispersions? 2 ^ We 
here follow Ref. [ll| and model the phonons using elastic 
continuum theory?^ 7 - Since even at room temperature, one 
effectively probes low energy scales, we keep only long- 
wavelength acoustic phonon modes. For these, previous 
work on related materials has show n 28 ' 29 that isotropic 
elastic continuum theory provides a reasonable approxi- 
mation. The phonon eigenmodes in the thin film geom- 
etry and their coupling to electronic modes have previ- 
ously been determined in the context of semiconductor 
quantum well structures^ (Note that the semi-infinite 
case has been treated in Ref. [H].) We basically repro- 
duce the phonon eigenmodes of Ref. but the coupling 
to the helical electronic eigenstates in a TI film is differ- 
ent from the semiconductor case. Note that piezoelectric 
couplings are suppressed by symmetry here?^ 9 - and spin- 
phonon type couplings 3 -^ are also expected to be subdom- 
inant to the deformation potential taken into account be- 
low. 

Most TI experiments have so far addressed only op- 
tical phonons? 3 -^ cf. also the corresponding situation 



2 



for Bi surfaces^ but very recently ARPES studies re- 
ported phonon-induced broadening of the lineshape in 
TIs.— ~— The observed Bi2Se3 electron-phonon coupling 
strength^ which has been extracted from the prefactor 
in the high-temperature quasiparticle decay rate T oc T, 
is in good agreement with our theoretical estimates.™ 
This indicates that the low-energy approach indeed pro- 
vides a reasonable starting point. To the best of our 
knowledge, no detailed measurements for the tempera- 
ture dependence of the TI film resistivity have been re- 
ported so far. We mention in passing that for the related 
case of a 2D graphene monolayer, a similar comparison of 
theor y 35 ' 36 to experimental has turned out to be success- 
ful. Remarkably, the electron-phonon coupling observed 
in Ref. and independently estimated by us^i turns out 
to be quite large. Under room temperature conditions, 
the resulting lifetime of helical quasiparticles is therefore 
short, and the resistivity is rather large. This behavior 
is substantially different from what is found in graphene. 
We suspect that this is (partially) due to the different 
Debye temperatures in both materials. 

The structure of the remainder of this paper is as fol- 
lows. In Sec. |TT] we discuss the model for the surface 
states in the thin film and their coupling to the quan- 
tized phonon modes. We then turn to the calculation of 
the electrical resistivity in Sec. IIII1 followed by the life- 
time broadening in Sec. IIV1 The paper concludes with a 
brief discussion in Sec. [V] Technical details of our cal- 
culations can be found in various appendices. Note that 
we use units with fi = 1 . 



II. MODEL 

In this section we describe the model employed in our 
study of electron-phonon scattering in a TI film. The 
model parameters below are chosen for Bi2Se3 as a con- 
crete example. The film has infinite extension in the xy 
plane and the width L, where \z\ < L/2. We start by re- 
viewing the construction of the effective surface Hamilto- 
nian describing the (upper and lower) electronic surface 
states of a TI film. 



A. Electronic surface states 

Keeping all terms up to second order in the momentum 
around the T point, (k x ,k y ,k z ), the low-energy physics 
of 3D TI materials like Bi2Sc3 or Bi2Tc3 is well de- 
scribed by an effective four-band modeli Using the ba- 
sis states {|P1+ t>> 1^2", t>, \Plt,±), 1^2",+)}, the low- 
energy bulk Hamiltonian reads^ — 







A,k. 



H 



<e + M -iA x k z 

iA x k z e - M A 2 k- 

A 2 k + e + M -iA 1 k z 

A 2 k+ iAik z e — M/ 



(1) 



with e =C + D^l + D 2 (k 2 x + k 2 y ), M = M - B x k\ - 
B 2 (k 2 + ky) and k± = k x ± ik y . The model parameters 
for Bi2Se3 have been determined from first principles,— 



M = 0.28 eV, C = -0.0083 eV, 

Ai = 2.26 eVA, A 2 = 3.33 eVA, 

B 1 = 6.86 eVA 2 , B 2 = 44.5 eVA 2 , 

Di = 5.74 eVA, D 2 = 30.4 eVA 2 . 



(2) 



We may write the Hamiltonian (fTJ) in the form H = Hq - 
%(k,) 
ho(k t ). 



if', where Hq 



is the 2x2 block 



matrix obtained for k x = k y = 0, with 



ho(k z ) = 



e {k z )+M Q -B 1 kl 
iAxk z 



-iAik z 



e (k z ) - M + B x k\ 



(3) 

Note that eigenstates of Ho have conserved spin. 

In order to find the surface states in the film geometry, 
we follow the usual strategy^— and first look for general 
bispinor eigenstates of ho, 

h (k z ^-id z )^(z)=E ^(z). (4) 

The general solution to Eq. (j4j) reads (j = ±, s = ±) 



*(*) = £ 



OS 



E a - C + M + (D 1 + BJrf 
-sAxrjj 



with arbitrary Cj s and the inverse lcngthscales 



(5) 



V± = 



—B ± Vb 2 - 4AC) /(2A) 



1/2 



where A = D\-B\,B = A?-2[M Bi+Di(C-Bo)], and 
C = (E - C) 2 - Mq 2 . The Dirichlet boundary conditions 
defining the film geometry, ^{z = ±L/2) = 0, then imply 
the transcendental equation 

[Eq - C + Mq + (D 1 + gi)?7 + ]?7- _ tanh(^L/2) 
[Eo - C + M + (£>i + Bi)rjl ]v+ ~ tanh(77+L/2) ' 

or the same condition with 77+ O r\- on the right hand 
side. Numerical solution of these equations yields the 
r point energies Eq '. The corresponding eigenstates 
*±(z) follow from Eq. (0, 



±[) J ^ ± { A 1 Fi(z) )' 

where the Af± are normalization constants and 
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(6) 



A. 
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77+ 00^(77+^/2) - i]- coth ± (?7_L/2) 



with coth + (y) = coth(y) and coth (y) = tanh(y). Fi- 
nally, the F functions are 



cosh (77+ 2) cosh(?7_z) 



cosh(77+L/2) cosh(77_L/2) 
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where F_ follows with cosh — ► sinh. Note that the eigen- 
states ^f±(z) describe both spin directions (a =t and 

(7=1). 

We now project the full Hamiltonian H to the ba- 
sis spanned by the surface states ([6]). We define Pauli 
matrices T a=x ^ y ^ z switching between the two solutions 
^t=±{z), Pauli matrices a a in spin space, and use To 
and do as identities. With the energy scales 

p(+) , p(-) 

the low-energy ("surface") Hamiltonian resulting from 
this projection reads 

H cS = E t <7 + ^T z ao - A 2 WT x (k x a x + k y a y ) + 0(k 2 ), 

(8) 

where W = ( 1 3/+|\t r _). The parameter A(L) is precisely 
the surface state gap described in the Introduction. For 
the parameters ©, -q _ is always real. However, 77+ is 
real for large L but purely imaginary for small L. In any 
case, we find that W is always real and positive. 

Noting that H e g commutes with r z cr z , it can readily 
be diagonalized by the unitary transformation U (k) = 
diag(U+, U-), where k = (k x ,k y ) and the U v= ±(k) are 
2x2 matrices in spin space, with v denoting the eigen- 
value of r z a z . With tana = 2A2W / |k|/A and tan# = 
kyjk x , we find 

(e' w / 2 cos(a/2) e~ w / 2 sin(a/2)\ 
Uv=+ ~ \-e ie ' % sm{a/2) e w ' 2 cos(a/2) ) ' W 

/-e- l9 / 2 sin(a/2) e^/ 2 cos(a/2)\ 
' ^ e i;e / 2 cos(a/2) e w l 2 sin(a/2) J' 

Switching to second-quantized notation, the eigenstatcs 
of -fTeff correspond to helical fcrmions with annihilation 
operator 

cu, vs = ^[U v (k.)]* as dk )T=tJ0 - i0 -, (10) 

(T 

where dk,rcr annihilates a spin-cr electron with in-plane 
momentum k in the transversal state \E' r (z). The low- 
energy electronic Hamiltonian (including the chemical 
potential fi) then takes the final form 

^el = ^ e k,s c ik,t;s C k,i;s: (H) 
k;ij,s— ± 

where the dispersion relation is 



ck,± = - £ °° - /" ± y a/1 + (2^2^/A) 2 k 2 . (12) 

We here choose the zero of energy by setting E§° = 

C + D 1 M /B 1 = liniL^oo E { ±] . For the parameters ©, 
we find Eq° ~ 0.22 cV. Moreover, for i — > oo, the length- 
scales ryjjj are given by T}7 ~ 12.3 A and ryZ 1 — 1.9 A. 
For fci ^> 1, the dispersion relation (fT2"|) is linear, with 
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Figure 1: (Color online) Electronic eigenstates for Bi2Se3 from 
Eqs. |T]) and ([2]). Main panel: Densities p T {z) in Eq. ([Til for 
L = 4 QL. Inset: Gap A vs thickness L. Note the semi- 
logarithmic scale. 

Fermi velocity vf — 2.77 x 10 5 m/s. Note that the index 
s = ± in Eq. pip does not correspond to spin anymore. 

Similarly, the particle density operator n(r, z) with r = 
{x,y) is written in terms of the dk^o- operators, 

h(r,z)= ^ r Pr{z)dl + ^ Ta d KTa . (13) 

k,q,r,CT 

Using Eq. (JTUJ), the density operator fTB")) can be trans- 
formed to the helical basis. We show the single-particle 
densities for the surface states [Eq. ||BJ], 

p T (z) = [#t . $ r ] ( z)) (14) 

in Fig. Q] for a film thickness of L = 4 QL, where 1 QL ~ 
9.5A for Bi2Se3.— This demonstrates that already for 
quite thin films, Eq. ^ describes surface states. Note 
that p r {z) is an even function of z. The inset of Fig. [1] 
shows the numerically obtained gap A(L), demonstrating 
the absence of oscillatory behavior for the parameters 
(J2J) as well as the exponential decay of A(L) due to the 
exponentially vanishing overlap of both surface states. 
We note in passing that for the parameters in Ref. [HI, 
Eq. instead predicts an oscillatory decay of A(L). 

B. Phonon model 

We now discuss the long-wavelength acoustic phonon 
modes in the TI film. Wc employ isotropic elastic con- 
tinuum theory, where the longitudinal (q) and trans- 
verse (c t ) sound velocities correspond to the two Lame 
constants.— In Bi 2 Se3, they are give n 39 ! 40 by q ~ 
2900 m/s and c t ~ 1700 m/s, respectively. Moreover, 
the mass density is pu = 7680 kg/m 3 .— In order to 
model the film geometry, we impose stress-free boundary 
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conditions^ at z = ±L/2. The quantized phonon eigen- 
modes for this problem have been determined in Ref. |3fj. 
For convenience, we briefly summarize the results next. 

Different phonon modes arc labeled by a set of quan- 
tum numbers, A = (q, A,n), where q = (q x ,q y ) is the 
surface momentum, A € (H, S, A) denotes the mode type, 
and n £ N is a branch index corresponding to the quan- 
tization of transverse momentum. The horizontal shear 
mode (A = H) decouples from all other modes and does 
not generate a deformation potential,— and we do not 
discuss this mode further. We are left with transversally 
symmetric (dilatational, A = S) and antisymmetric (flex- 
ural, A = A) phonons. Denoting the dispersion relation 
of a given phonon mode A by J7a (see below) and the 
surface area by A, the displacement field operator is 



U(r,M)=£ 



e i(q-r-Sl A t) 

^2 PM Att A 



u A (*) &a + H.c, (15) 



where 6a is a bosonic annihilation operator and the non- 
interacting phonon Hamiltonian is 



ff ph = fi A ( b A b A + 1/2 



(16) 



The orthonormal eigenmodes u\(z) in Eq. (|15|) describe 
linear combinations of e ±tk '-t z waves, where 



ki, t = ^(Oy/cm) 2 - q 2 ; 



(17) 



h,t = ini,t with Kij = \Jq 2 - (^a/q,*) 2 for Ha < a.tq- 
Writing ua(z) in the form 



u(z) 



iq<Pi - -j- J e q + ( + iqcj) t | 



d<pi 
dz 



(18) 



where e q = q/g and 

4>i,t = aij cos(fc; jt z) + bij sm(k^ t z), (19) 

the stress-free boundary conditions at z = ±L/2 yield 

d(f>i 
dz 



0. 



(20) 



2ig^ + ( g a-^ I = 0. 
dz 



Since both equations have to be fulfilled at z = ±L/2, we 
have four linear equations for the four unknown parame- 
ters {ai_ tl bi t ). Setting the corresponding determinant to 
zero, we obtain the following two possibilities. First, for 
symmetric modes (A = S), we have the condition 



(q 2 - k 2 ) 2 cos(fc ; i/2) sin(fc t i/2) 
Aq 2 hk t sin(fc ( i/2) cos{k t L/2) 



(21) 



0. 



Numerical solution of this transcendental equation ob- 
tains the quantized set of dilatational phonon frequen- 
cies ^A=(q,s,n)- The corresponding eigenvector, u\(z), 
follows from Eqs. ([T8"]) and (fT9"j) with a t = bi = and 



ai 



2Af S q 
cos(fe;i/2)' 



iM s (q 2 



k t cos(k t L/2) 



(22) 



& 

< 




Figure 2: (Color online) Phonon dispersion relation, J7a vs q, 
for the symmetric (A = S) mode (red solid curves). Shown 
are the ten lowest branches corresponding to the index n. 
Dashed lines separate regions I, II, and III (see main text). 
The dash-dotted line gives the dispersion relation in Eq. (|24|) ; 
note that the n — 1 mode coincides with the Rayleigh mode 
for qL > 1. 



Second, for antisymmetric modes (A = A), we arrive 
again at the condition in Eq. (f2~T|) but with the ex- 
change cos <-> sin. Solving that equation yields the 
set ^A=(q, J 4,n) °f quantized flcxural phonon modes. The 
eigenvector ua(z) follows again from Eqs. (fT5|) and (|19l) . 
where now ai — b t = and 



2A/kg 
sin(fc;i/2) : 



a i 



iM A (q 2 



kt sin(fc t L/2) 



(23) 



The normalization factors Af\=s.A appearing in Eqs. (|22|) 
and (|23| are given in Appendix [Al 

Numerical solution of Eq. (|2"Tj) yields the spectrum, 17a, 
for the symmetric mode (A = S). The result is shown in 
Fig. [21 We distinguish three different regions, namely a 
case where both ki and kt are purely imaginary (region I), 
when only ki is purely imaginary but kt is real (region II) , 
and finally a case where both fc; and k t are real (region 
III). We observe from Fig. [5] that the n = 1 mode is the 
finite-width analogue of the well-known Rayleigh surface 
modej 27 i 31 For the semi-infinite geometry, the Rayleigh 
mode is the lowest-lying phonon42 It has the dispersion 
relation 



il = c R q, c R ~ 0.92c t . 



(24) 



In fact, for qL ^> 1, both Eq. (f2"Tj) and the corresponding 
equation for A = A reduce to 

(q 2 + k 2 ) 2 = Aq 2 K t Ki. 

As discussed in Ref. H3, this equation readily yields the 
sound velocity cr of the Rayleigh mode. 
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C. Electron-phonon coupling 

The dominant coupling of the above phonon modes 
to the electronic surface states comes from the deforma- 
tion potential^ which couples the local electronic den- 
sity n(r, z) [Eq. (|T5|) ] to the divergence of the displace- 
ment vector, V • U(r, z), see Eq. (fTS")) . Since the surface 
state density p T {z) in Eq. (|14|) is even in z, the anti- 
symmetric phonon mode (A = A) does not couple to the 
surface states. We therefore keep only the symmetric 



phonon mode from now on (and omit the index A = S). 
Transforming Eq. (fTB")) to the helical basis, see Eq. (|10|) . 
the second-quantized electron-phonon coupling Hamilto- 
nian reads 

a s-^, ,(«,»,«'), f tt 

tle-ph — -7= 2^ k,q,« °q^ lC k+q,t;s C k,t;s' + W- -: 

q,k, n;i;,s,s' 

(25) 

where the M matrix elements involve the unitary matri- 
ces [U v (k)] sa in Eq. ©, 



J 



M, 



(^,s,s ; ) 



k,q,n 



ft,, 



^[^(k + q^C/^k)], 



L/2 



dz p T=va (z)(f>i(z), 



(26) 



-L/2 



r 



with the phonon dispersion ft g , n in Fig. [2] 0; is given by 
Eqs. P^|) and (JHJ). The deformation potential strength 
a in Eq. (|2"5]) can be estimated as follows. The high- 
temperature behavior of the on-shcll imaginary part of 
the electronic self-energy is (see Sec. IIV[) 



ImS(fc, T) = —TrXkksT, 



(27) 



which allows to experimentally extract the dimension- 
less effective electron-phonon coupling constant A&. The 
relation (f27f has been observed for Bi2Sc3 in ARPES 
experiments^ and A = 0.25±0.05 has been measured. In 
these experiments, the Fermi level was near the bottom 
of the conduction band, p ~ 0.28 eV, and k in Eq. (|27|) 
corresponds to energies ~ 50 to 100 meV above the Dirac 
point. Computing Xk within our model, see Sec. [TV] the 
observed value for A corresponds to a = (30 ± 8) eV. We 
employ the value a = 30 cV below. 

The total Hamiltonian employed in the following sec- 
tions is then given by H = H c \ + i/ p h + _ff c _ p h, see 
Eqs. ([TT]), (H6l) and ([2%]). We first address the phonon- 
induced resistivity p in Sec. IIIII and then turn to the 
quasiparticlc lifetime in Sec. IIVI 



the Hamiltonian described in Sec. HU As explained in 
Sec. Ill C[ only symmetric (dilatational) phonon modes 
can cause a finite resistivity for the low-energy sur- 
face states within the bulk gap. We compute p within 
the framework of the linearized Boltzmann equation^ 
which has also been employed previously for the related 
graphene case. 35 ' 36 The resulting quasiclassical estimate 
for p is valid 3 - 5 - as long as p is small compared to the re- 
sistance quantum, p <C h/e 2 ~ 25.8 kO. We sketch the 
standard derivatio n 30 ' 35 ' 36 ' 44 for p in Appendix [B] The 
result takes the form 



1 

P 



- y 
2 ^ 

v,s—± 



dk 



vl s T v (e k>s )[-d e n F (ek,s)}, (28) 



III. RESISTIVITY 

Here we discuss the T-dependent phonon contribu- 
tion to the electrical resistivity, p, in the TI film, using 



where the dispersion relation for helical fcrmions 
[Eq. (TT21 ] defines the group velocity, v^ s — ik ■ Vk£k,s- 
Moreover, ni?(e) is the Fermi function, and the energy- 
dependent electron-phonon transport scattering rate (in- 
verse time) is 



J 



q,s' 



1 - 



Vk+q,s> 



COS( 



k,q 



1 - n F (e k+q;S /) 
1 - n F {ek tS ) 



(29) 
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where #k.q is the angle between k and k + q, and the transition probabilities are obtained from Fermi's golden rule. 
Using Eq. ([2"5"|). we find 

2 

5 (eic,a + Vftq,n ~ £k+q,s') , (30) 

i 



ir. 



(k,^s)— >(k+q,ws') 



E 



M, 



(s',s) 
k,q,n 



where 715(e) is the Bose function. While the M matrix 
elements (|26p depend on the index u = ±, we note that 
\M\ 2 and therefore the transition probabilities W are v- 
independent. This also implies that t v does actually not 
depend on v. 

With the polar angle between k and q, such that 

k + qcosd 

COSflkq = -, (31 

y/k 2 + q 2 +2kqcos9 
I 



the angular integration in Eq. (j2"9")l can be encapsulated 
in the "transport Eliashberg function" , see also Ref. [TJ 



Q " L 



(V, s ) 

k.q.n 



(Ck,s + Z^q.ra — £k+q,s') 



This allows us to write the momentum relaxation rate (|29[) in the form 

1 



~0k, s ) 



/ *^;, s (9)^B(^,,n) : ~ , i 

Jn 1 — "FlCksj 



(32) 



(33) 



The ^-integration in Eq. (j3"2")l can then be carried out an- 
alytically. We quote the (lengthy) result in Appendix ICl 
which is useful when computing T numerically. For low 
temperatures, the quasi-elastic approximation, £l q _ n <C 
y/ (A/2) 2 + (A2Wk) 2 , is applicable and allows to sim- 
plify the full result for T to the v- independent form 



•Ffe.n.sC?) = 6(2fc-g) 



1 



(34) 



V(A/2) 2 + (A 2 Wkf 
(A 2 Wk) 2 



M, 



k,q,n 



where 8 = 6q (see Appendix [Cj determines the polar 
angle between k and q appearing in the matrix element 
M, and the Heaviside function is denoted by @(y). Note 
that there is no contribution from interband transitions 
at low temperatures. 

The crossover temperature from the low- to the high- 
temperature behavior in this system is se t 11 ' 35 by the 
Bloch-Griincisen temperature, 



Tbg = 2k F c R /k B , 



(35) 



with the Rayleigh velocity cr in Eq. (|24|) . kp(L) is de- 
fined by ek F ,s=+ = with the dispersion relation (|T2l . 
For T <C TrGi the J- function can be approximated 
by the quasi-elastic expression [Eq. It receives 



the dominant contribution from the n = 1 branch cor- 
responding to the Rayleigh surface phonon. For small 
q, we find f2 9l n=i = c s q with c s = 2754 m/s (which is 
slightly below q), see also Fig. [21 In addition, we have 
(pi(z) = 2(c t /c s ) 2 /(q\/Z) and c s < min(|u fcji ,|, A 2 W), 
leading to 



ft 



(9) 



(c t /Q) 4 



This allows us to perform all remaining integrations and 
yields a T 4 law for the resistivity at low temperatures, 



A 



p(T « T BG ) = -^A 



8jkpa 2 



T 

e" \ Tbg 
where the dimensionless prefactor A is 

ClCg J V 
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np M v% F c s 
rf.r 



(36) 



(37) 



2e x 



[(n 2 +x 2 ){e x + l)] 2 



68.4295. 



For L — > oo, A obviously vanishes. This suggests that 
for elevated temperatures (but still T < T BG ) and finite 
L, the T 4 law is replaced by the i-independent p cx T 5 
law found in Ref. [Ill We can estimate the crossover 
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Figure 3: (Color online) Phonon contribution to the resistiv- 
ity p vs temperature T for a TI film of width L — 4 QL and 
several values of the chemical potential p. Dashed lines indi- 
cate the analytical results for low [Eq. (|36l) ] and high [Eq. (|38l) ] 
temperatures. Note the double-logarithmic scale. 




5 10 15 20 

L (QL) 

Figure 4: (Color online) Width (L) dependence of the phonon 
contribution to the resistivity p for p = 0.2 eV and sev- 
eral temperatures. The dashed horizontal line indicates one- 
quarter of the resistivity poo(T) in the semi-infinite geometry 
with otherwise identical parameters.— 



temperature T c as follows. For T < Tbg, we expect an 
expansion of the form 

(e 2 /h)p = A(T/T BG f + |(T/T BG ) 5 , 

with A cx 1/L in Eq. (|36l) and the L-independent con- 
stant B given in Ref. [ill The crossover from the T 4 law 
(for T < T c ) to the T 5 law (for T c < T < T BG ) thus hap- 
pens around the temperature T c = (4A/B)Tbg- This 
gives T c ~ 0.14TBG/(fc_F^), which is independent of the 
chemical potential since Tbg kp. For L = 4 QL, we 
obtain T c w 0.9 K. The T 4 law can thus only be observed 
for very thin and clean TI films. 

In the opposite high-temperature limit, essentially all 
phonon branches indexed by n contribute to the trans- 
port Eliashbcrg function (pl2"j) . see Appendix C. Then the 
relaxation rate T _1 (ek, s ) in Eq. (|3"3]) is basically a linear 
function of the energy. Since the linear term does not 
contribute to p after integration in Eq. (|28p . we obtain 
the approximation 1/p ~ (e 2 /h)vk F kFT(e = 0), where 
Eq. (|33|) yields the linear high-temperature law 



(38) 



p(T » Tbg) = \c^- 
e. Jbg 

with the dimensionless prefactor 



c=^- E r qdq 



(39) 



q.n 



Next wc show the full temperature dependence of p 
obtained numerically for a fixed width L = 4 QL and 
several values of the chemical potential p, see Fig. [3] 



In that case, when measured relative to E§°, we have 
Eq ~ 16 meV and A/2 ~ 13 meV. For the lowest ^ in 
Fig. [3l the Fermi level is thus located inside the surface 
gap and one has a very large resistivity, where the quasi- 
classical approach is not reliable in any case. For low 
temperatures, T < Tbg, the analytical result ([55)1 with 
p oc T 4 is nicely reproduced by numerics. In this temper- 
ature regime, only the Rayleigh mode (n = 1) is relevant, 
similar to what one finds in the semi-infinite geometryi^ 1 - 
In the high-temperature limit, both the p oc T scaling 
and the prefactor C in Eq. (f3"9")) are also consistent with 
our numerical findings. 

Finally, Fig. @] shows the width (L) dependence of p 
at fixed chemical potential and for several T. Two note- 
worthy observations can be drawn from Fig. 01 First, for 
low temperatures we observe a "dip" in Fig. 21 where 
p{L) < p(L — > oo ) for intermediate values of L. Second, 
for L — >• oo, p(L) approaches 1/4 of the single-surface 
value Poo{T) obtained for the semi-infinite geometry* 1 ^ 
Naively, we would expect p(L — > oo) = Poo/2 because of 
the presence of two surfaces in the film geometry. This 
discrepancy indicates that the L — > oo limit is singular, 
and it is not possible to really decouple both surfaces in 
such an interacting system; see also Ref. [2l| for a related 
discussion. 



IV. LIFETIME BROADENING 

Next we discuss the quasiparticle lifetime (inverse de- 
cay rate) for the surface fermions in the TI film due 
to their coupling to phonons, see -ff e - P h in Eq. (pj"5"j) . 
which implies a finite linewidth of ARPES spectral fea- 
tures. The decay rate, Tk(T) = —2 Im S, follows from the 
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"rainbow" diagram yields the self-energy 



Figure 5: (Color online) Main panel: T-dependence of the 
decay rate T of a TI film of width L = 4 QL for k = k F and 
k — 0.5k F . For k — 0.5k F , only the \i — 0.2 eV result is dis- 
played. Dashed lines indicate the low- and high-temperature 
laws (F oc T 2 and oc T), respectively. Inset: fc-dependence of 
r for \x = 0.2 eV and two different temperatures: T = 3 K 
(solid line) and T = 300 K (dashed line; the shown result has 
to be multiplied by 10). 



0.48 



0.46 



0.44 



0.42 



0.40 



Figure 6: (Color online) Width (L) dependence of the effective 
electron-phonon coupling constant at the Fermi level Afc F for 
fi = 0.2 eV. The dashed horizontal line indicates one-half of 
the effective coupling constant in the semi-infinite geometry 
with otherwise identical parameters.— 




dq 

f 



(s',s) 
k,q,n 



(40) 



E 

v=± 



LU + i.Q+ + lsQc^n - Ek+q.s' 



Introducing the Eliashberg function (q) exactly as 

the transport Eliashberg function T in Eq. ([52)1 but with- 
out the factor [1 — (fk+q,s'/ u k,s) cos#k,q]j the quasipar- 
ticle decay rate follows as 



r*(T) 



E 



(41) 



x [n B {Q q ,n) + n F (Vt q ^ n + i>£fc,+)]. 



Expanding this result for high temperatures, T Tbg, 
as in Sec. |TTT] yields, see also Eq. (f2"Tf. a linear In- 
dependence, 



r fc (r»T BG ) = 2nX k k B T, 



(42) 



A fc = — 



2tt ^ 



qdq- 



F, 



+ («) 



The L-dependence of A& is shown for k = k F in Fig. |6] 
We observe an oscillatory dependence, with a saturation 
at one-half of the corresponding semi- infinite result. 

For low temperatures and k = k F , the decay rate is 
dominated by the n = 1 phonon mode with q — > 0. After 
some algebra, we find that this implies a T 2 law, 



r fcF (T«T BG ) 



Mct/ci) 4 (k F c R a) 2 1 



T 



(43) 

Again, when T > T c , the T 2 law (which scales oc 
1/L ) co mpetes with the i-independent T 3 law found in 
Ref.TTll. see Sec. [Til] Finally, when k ^ k F and T < T BG , 
the quasiparticle decay rate saturates at the finite value 



/■OO 

IVfc F =« 2 E / 9 d 9 (l £ fc+l 
„ Jo 



with = sgn(fcjr — fc). 

Figure [HI shows that the £ — > oo limit of the decay rate 
always tends to r oo (T)/2, where is the corresponding 
decay rate for the semi-infinite geometry^ This discrep- 
ancy with the naive expectation T(L — > oo) = has the 
same origin as the anomalous factor 1/2 appearing in the 
large- L behavior of the resistivity discussed in Sec. Mil 



V. CONCLUSIONS 



imaginary part of the on-shell self-energy E s=+ (k, tu = 
e k.s=+)^ Expanding up to second order in _ff _ph, the 



In this paper we have studied the effects of long- 
wavelength acoustic phonons on the topologically pro- 
tected surface fermions in topological insulator films. 
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Our model employs the established low-energy electronic 
Hamiltonian and an isotropic elastic continuum approach 
for the phonons, with the deformation coupling provid- 
ing the dominant interaction mechanism. The electron- 
phonon coupling turns out to be surprisingly strong, in 
accordance with recent ARPES results.— 

Using a quasiclassical approach, we have computed 
the temperature-dependent resistivity of the film due to 
phonon backscattering, and found a linear T dependence 
above the Bloch-Griineisen temperature. In this temper- 
ature regime, the phonon-induced resistivity can over- 
come the disorder-induced (T-independent) contribution 
and should be observable with present samples. Simi- 
larly, the linear T dependence of the quasiparticle decay 
rate found here is observable 2 ^ in ARPES experiments. 
The low-temperature behaviors of the resistivity and of 
the quasiparticle decay rate are probably more difficult 
to observe. 

An interesting extension of our work would be to in- 
clude the effects of a magnetic field. Magnetotransport 
measurements in thin films were recently performed^ 6 - 



and found clear evidence for Landau level formation asso- 
ciated with the massless Dirac fermions forming on both 
surfaces. The observed broadening of the Landau lev- 
els was assigned to disorder and/or interaction effects, 
but at elevated temperatures, our analysis indicates that 
electron-phonon interactions may be relevant as well. 
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Appendix A: Normalization constants 

Here we provide the normalization constants Ms, a ap- 
pearing in Eqs. (JUJ) and {23) in Sec. Ill Bl Specifically, 
we get these constants after some algebra from 



Ms 2 



2c? 
2c 2 



4q 2 L 



_cf cos 2 (fc/L/2) 



>2 



4q 2 L 



cf sin 2 (fc ; L/2) 



sin(fc;L) 

sin(fc;L) 
k,L 



(q 2 ~ k?) 2 L ( S m(k t L) 

k 2 cos 2 (k t L/2) V hL 

(q 2 -k 2 ) 2 L 
k 2 sm 2 (k t L/2) 



sm(k t L) 
k t L 



-8( g 2 -fc t 2 ) tan(W2) " 
h 

+ 8(g 2 -fc 2 ) C ° t(fctL/2) 



hi 



Appendix B: Linearized Boltzmann equation 

The linearized Boltzmann equation has been derived 
for closely related problems before , 30 ' 35 ' 36 ' 44 and we here 
follow those works and briefly sketch the derivation of 
Eq. (|28[) . In the quasiclassical approximation, the quasi- 
particle distribution function /(r,k, t) (for simplicity, we 
here omit the v and s indices) obeys the well-known 
Boltzmann equation. In the absence of a force F = — eE 
due to the external electric field, / reduces to a Fermi 
function, / = n^(e k ), and the collision integral vanishes. 
In the presence of the force, / is expanded in terms of 
Lcgcndre polynomials P„(cos a), where a is the angle be- 
tween k and F. Keeping only terms linear in F, we have 
/(k) = ni?(ek) + cos(a)/i(e k ). Using detailed balance, 
for given transition matrix elements Wk-^k', we obtain 
the linearized Boltzmann equation (LBE), 



F • v k 9 e ?i F (e k ) = 
np(ek) 



E 

k' 



n F {ek>) 



cos(o/)/i(ek') 



WW' 

1 - n F (ew) 



1 - n F {ek) 



cos(a)/i(e k ) 



tion 
1 



r(e k ) 



k' 



VUk- 



^k' T(e k ') 
v-k r(e k ) 



cos(z9) 



1 - n F (ey) 
1 - n F (e k ) ' 



where i? is the angle between k and k'. The solu- 
tion for r(e k ) determines the electron momentum re- 
laxation time. When the scattering of quasiparticlcs 
from long-wavelength acoustic phonons is quasielastic, 
Qq,n "C \/j.\, we can set T(e k ') = T(e k ) for the right-hand 
side of the above integral equation; this is equivalent to 
the "test particle approximation."— The current density, 
j = — (e/A) Y^ik v k/(k), points parallel to the electric 
field direction and has the magnitude 

e 2 E 

j = — ^w k cos 2 (a)r(e k )[-a £ n f (e k )]. 
* k 

Using Vk = Vk and performing the angular integration, 
we arrive at the phonon contribution to the resistivity 
quoted in Eqs. {25]) and {29]). 



Appendix C: Transport Eliashberg function 



Using the Ansatz /i(e k ) = — r(e] SL )v] SL Fd e n F (e k ), after 
some algebra the LBE leads to the linear integral equa- 



We here give the analytical result for the full trans- 
port Eliashberg function T defined in Eq. (|52^) . Some 
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straightforward yet tedious algebra allows to perform the ^-integration. We find the (lengthy) result 
I 



A 



» 



e 



k-q)Q(q 



Q 



k,q,n,s 



ir(A 2 Wy 



{i 2 -(Q ( kLs-kr) ((Qi:Us + kr- q 2 ) 



1/2 



e 



.4 



k,q,n,s 



A/2 



k,q,n,s 



l-i 



2 g 2 



k,q,n,s , 



A 



2fc 2 



k.q.n 



where we use the notations 
iM 



.4 



k,q,n,s 



^ ks,q,n,s 



V(A/2) 2 + {A 2 Wkf + vn q , n , 



AoW 



and the polar angle 9 = 6 e [0,7r] follows from 

y/h?+q* + 2kqCO&0 Q = QS,„ jS) 

fixing the polar angle between k and q in the matrix 
element M. 
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